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Abstract 

The past decade has seen the rise of various radio astronomy arrays, particularly for low-frequency 
observations below lOOMHz. These developments have been primarily driven by interesting and fun¬ 
damental scientific questions, such as studying the dark ages and epoch of re-ionization, by detecting 
the highly red-shifted 21cm line emission. However, Earth-based radio astronomy below frequencies 
of 30MHz is severely restricted due to man-made interference, ionospheric distortion and almost 
complete non-transparency of the ionosphere below lOMHz. Therefore, this narrow spectral band 
remains possibly the last unexplored frequency range in radio astronomy. A straightforward solution 
to study the universe at these frequencies is to deploy a space-based antenna array far away from 
Earths’ ionosphere. Various studies in the past were principally limited by technology and computing 
resources, however current processing and communication trends indicate otherwise. Eurthermore, 
successful missions which mapped the sky in this frequency regime, such as the lunar orbiter RAE-2, 
were restricted by very poor spatial resolution. Recently concluded studies, such as DARIS (Dis- 
turbuted Aperture Array for Radio Astronomy In Space) have shown the ready feasibility of a 
9—satellite constellation using off the shelf components. The aim of this article is to discuss the cur¬ 
rent trends and technologies towards the feasibility of a space-based aperture array for astronomical 
observations in the Ultra-Long Wavelength (ULW) regime of greater than 10m i.e., below 30MHz. 
We briefly present the achievable science cases, and discuss the system design for selected scenarios, 
such as extra-galactic surveys. An extensive discussion is presented on various sub-systems of the 
potential satellite array, such as radio astronomical antenna design, the on-board signal processing, 
communication architectures and joint space-time estimation of the satellite network. In light of 
a scalable array and to avert single point of failure, we propose both centralized and distributed 
solutions for the ULW space-based array. We highlight the benefits of various deployment locations 
and summarize the technological challenges for future space-based radio arrays. 
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1 Introduction 


The success of Earth-based radio astronomy in the frequencies between 30MHz and 3GHz, is jointly 
credited to Earth’s transparent ionosphere and the steady technological advancements during the past 
few decades. In recent times, radio astronomy has seen the advent of a large suit of radio telescopes, 
particularly towards the longer observational wavelengths i.e., > 3m. These arrays include the Murchi¬ 
son Widefield Array (MWA) (Lonsdale et ah, 2009), LOw Erequency Array (LOEAR) (van Haarlem 
et ah, 2013) and the Long Wavelength Array (LWA) (Ellingson et ah, 2009) to name a few. However, 
Earth-based astronomical observations at these ultra-long wavelengths are severely restricted (Kaiser and 
Weiler, 2000). Eirstly, due to ionospheric distortion, especially during the solar maximum period, when 
scintillation occurs and the celestial signals suffer from de-correlation among the elements of a ground 
based telescope array (Kassim et ah, 1993). Advanced calibration and mitigation techniques which are 
currently employed in LOEAR telescope array, can be used to remove these distortions, provided the 
time scales of disturbances is much longer than the time needed for calibration process (Wijnholds et ah, 
2010). Eurthermore, at frequencies below lOMHz the ionosphere is completely non-transparent which 
impede observations by ground-based instruments. In addition to ionospheric interference, man-made 
transmitter signals below 30MHz also impede astronomical observations. This terrestrial interference was 
even observed as far as ~400, 000km away from Earth by the RAE-2 lunar orbiter, which was limited by 
very poor spatial resolution at these wavelengths, e.g., 37° at 9.18MHz (Alexander et ah, 1974). Due 
to the above mentioned reasons, the very low frequency range of 0.3 — 30 MHz remains one of the last 
unexplored frontier in astronomy. An unequivocal solution to observe the radio sky at ULW with the de¬ 
sired resolution and sensitivity is to deploy a dedicated satellite array in outer-space. Such a space-based 
array must be deployed sufficiently far away from Earths’ ionosphere, to avoid terrestrial-interference 
and offer stable conditions for calibration during scientific observations. 

1.1 Science at ultra-long wavelengths 

A space-based low frequency radio instrument would open up the virtually unexplored ULW domain 
and as such addresses a wealth of science cases that undoubtedly will lead to new exciting scientific 
discoveries, similar to the uncovering of other wavelength domains has revealed in the past. Eor most of 
these science cases such an array would add information to the existing radio, optical, infrared, sub-mm 
or high frequency X-ray or gamma-ray instrumentation, and thereby providing insight into the processes 
that take place at the lowest energies and largest physical scales. As explained in (Jester and Ealcke, 
2009) and (Klein-Wolt et ah, 2012) science topics include for example the study of the solar activity and 
space weather providing important clues on the effect of solar flares and bursts on the Earth to much 
larger distances from the solar surface. Investigating the magnetospheric emission from large planets 
such as Jupiter and Saturn (see Zarka et ah, 2012) reveals information on the spin period of these 
planets. Other key science cases include the study of large-scale structures from galaxy clusters (and 
radio galaxies) and the detection of Jupiter-like flares and Crab-like pulses from (extra-)galactic sources. 

The greatest advance in science is expected in the study of the very early universe in a period referred 
to as the cosmological Dark Ages (Rees, 1999). The only direct window to this period in the evolution of 
the early universe is provided by the 21-cm line from the neutral hydrogen caused by the spin-flip of the 
electron which, redshifted up to a factor 1000, is now visible in the low frequency radio regime between 
1.4 - 140 MHz. The Dark Ages is the period between the epoch of recombination when the universe 
became transparent and the cosmic (microwave) background was emitted, and the epoch of reionisation 
(EoR) when the flrst stars started to reionise the neutral hydrogen. The hydrogen has played a major 
role in the Dark Ages and the EoR, being the only source of light in the absence of stars and through 
couplings with the cosmic background radiation and the surrounding neutral hydrogen it allows us to 
trace the distribution of matter in the early universe and observe the formation of the very flrst stars 
and structures in the universe (Ciardi and Eerrara, 2005). The global Dark Ages signal, essentially the 
redshifted 21-cm line absorption feature, is expected to peak around 30-40 MHz and is weak, ^10^ below 
the foreground signal. 

However, (Jester and Ealcke, 2009) show that with a single antenna placed at an ideal location on the 
moon (i.e. under low REI and stable temperature and gain conditions) the global signal can be detected 
at a 5 cr level for one year of integration. In order to trace the variations in the hydrogen at arc-minute 
or even arc second scale resolution a larger sensitivity and hence collecting area is required. Eor instance. 
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(Jester and Falcke, 2009) show that up to 105 individual antenna elements are required corresponding 
to 0.5 km^ in order to reach 10' spatial resolution (see Loeb and Zaldarriaga, 2004). So, in addition to a 
low-RFI, low-temperature and stable gain location, a ULW radio interferometer that aims at detecting 
now only the global Dark Ages signal but also the arcmin variations in the mass distribution of the Dark 
Ages and the EoR must have a large collecting area (in the order of km^). Ideal locations include the 
farside of the moon, an eternally dark crater on the lunar south or north pole, or space-based solutions 
such as for instance in a Sun-leading or trailing orbit, at the Sun-Earth L2 point or in Lunar orbit. 

1.2 Previous studies 

The proposition for a space-based radio astronomy instrument is not novel (Weiler et ah, 1988; Basart 
et al., 1997a,b; Kaiser and Weiler, 2000). One of the first such proposal was by Gorgolewski (1965), 
who discussed the benefits of a moon-based radio interferometer. In 1973 the lunar orbiter RAE-2 was 
launched, which mapped the non-thermal galactic emission in the frequency range of 25kHz to 13MHz 
using a 37m dipole antenna and achieving a resolution of 37° (Alexander et ah, 1974). Science at the 
Ultra-long wavelengths was revived from 1980s with a particular focus on Lunar based arrays (Burke, 
1990; Burns et al., 1990). The Lunar surface on the far-side presents a large and stable platform for 
antennas and shields unwanted interference from Earth and the Sun (Woan, 1999; Kuiper and Jones, 
2000; Takahashi, 2003; Aminaei et ah, 2014), which motivated studies such as VLEA (Smith, 1990), 
MERIT (Jones et ah, 2007) and more recently DEX (Klein-Wolt M. et ah, 2013). Along similar lines, 
lunar orbiting single-satellite missions dedicated for radio astronomy such as LORAE (Burns, 1990) and 
DARE (Burns et ah, 2012) were also investigated to map bright sources and to facilitate relatively easier 
Earth-based down-link of science data. Eurthermore, the pursuit of higher angular resolutions has led 
to Earth-orbiting single-satellite missions such as HALCA (Hirabayashi et ah, 2000) and Radio Astron 
(Kardashev et al., 2013) which enable Earth-space very long baseline interferometry. 

The concept of space based ULW array for radio astronomy however has seldom been explored 
adequately, which is our focus in this article. A notable study in this regard was the ALEA concept, 
which proposed an array of 3 — 16 satellites in the radio quiet L2 points as a deployment location (Jones 
et al., 1998). More recently, two ESA funded studies namely EIRST (Bergman et ah, 2009) and DARIS 
(Boonstra et ah, 2010) investigated passive-formation flying missions for space-based satellite arrays (see 
Table 1). The EIRST study proposed a constellation of 7-satelhtes deployed at the second Earth-Moon 
Lagrange (L2) point, which allowed for a low-drift orbit and yet remained sufficiently far enough from 
Earth to avert interference. On the other hand, the DARIS mission primarily investigated the feasible 
ULW science cases and showed ready feasibility of 9-satellites using existing off the shelf technologies. 
The benefits of both the studies are combined in the SURO concept, which proposes a mission with 
a larger observational frequency range at Sun-Earth L2. In all these studies, a dedicated centralized 
mothership managed the processing and communication. However, futuristic arrays such as OLEAR 
with >10 satellites will operate cooperatively and employ a distributed architecture for both processing 
and communication. 

1.3 Overview 

The purpose of this article is to discuss the current technological advances towards the feasibility of space- 
based array for radio astronomy at ultra-long wavelengths. To this end, we elaborate on the system design 
for a space-based array in Section 2. We address various subsystems of the potential satellite array in the 
Sections 1.1 - 6, including the astronomy antenna design in Section 3. While current technologies limits 
us to < 10 nodes, we foresee next generation arrays will contain larger number of satellites and operate 
as a co-operative wireless network. Hence, a dominant theme of the article is to investigate the extension 
of the proposed centralized solutions to distributed scenarios, particularly for processing (Section 4), 
communication (Section 5) and joint space-time estimation of the satellites in the network (Section 6). 
We summarize the article with a brief overview of the potential deployment locations (Section 7) and 
the fundamental challenges ahead for a space-based ULW array (Section 8). 
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FIRST 

DARIS 

SURO-LC 

OLFAR 

Timeline 

2009 - 2010 

2009 - 2010 

2011 - 2012 

2010 - 2014 

No. of satellites {N) 

6t+lt 

8t + lt 

> 8t + it 

> 10 

No. of polarizations (Np^i) 

3 

3 

3 

3 

Obs. frequency (ly) 

0.3 - 50 MHz 

0.3 - 10 MHz 

0.5 - 60 MHz 

0.3 - 30 MHz 

Instantaneous BW (Aza) 

lOOKHz 

1 MHz 

IMHz 

> 1 MHz 

Obs. wavelength (A) 

600 — 6 m 

10^ - 30 m 

600 — 5m 

10^ - 10m 

Longest baseline 

30 km 

100 km 

30 km 

100 km 

Spatial resolution 

0.6 at 50MHz 

1' at lOMHz 

0.5' at 60MHz 

0.3' at 30MHz 

Array architecture 

Centralized 

Centralized 

Centralized 

Distributed 

Estimated Mass 

- 200Kgt, lOKg^ 

- 550Kgt, lOOKg^ 

- 500Kgt, lOKgt 

< 5Kg 

Deployment location(s) 

Sun-Earth L2 

Dynamic solar orbit. 
Moon far side, 
Sun-Earth L2 

Sun-Earth L2 

Earth/Moon orbit 


Table 1: Recent Space-based aperture array studies: An overview of system requirements for various space-based 
aperture array feasibility studies for ultra-long wavelength observations, namely DARIS (Saks et ah, 2010; Boonstra et ah, 
2010), SURO-LC (Baan, 2012) and OLFAR (Bentum et ah, 2009; Rajan et ah, 2011), where f and J denote mothership 
and daughter node respectively. 


2 Ultra-long wavelength interferometry 

2.1 Aperture synthesis 

Radio astronomy imaging is achieved by aperture synthesis, where in the cosmic signals received at a 
large number of time-varying antenna positions, are coherently combined to produce high quality sky 
maps. For a A/"—antenna array, each antenna pair forms a baseline of an aperture synthesis interferometer, 
contributing N = 0.bN{N — 1) unique sampling points at a given time instant. Let Xi(t) and Xj(t) be 
two arbitrary antenna position vectors at time t forming a baseline, then the corresponding uvw point 
is defined as 

[uij{t),Vij{t),Wij{t)\= (Xi(i)-Xj(i))/A, (1) 

where A is the observed wavelength. Figure 1(a) shows the uvw for a A/" = 9 satellite cluster which is 
arbitrarily deployed with a maximum distance separation of d = 50km and an observational frequency 
of lOMHz. The effective synthesized aperture is then obtained by projecting the uvw points onto a 
2-D plane which is orthogonal to the source direction. As an illustration. Figure 1(a) shows 3 such 
projections (in black) for sources orthogonal to the uv^ uw and wv planes. The minimum distance 
between the satellites is only constrained by practical safety requirements and the maximum distance 
between the satellites defines the resolution of the interferometric array as 

0 = \/d. (2) 

The Van Cittert Zernike theorem relates the spatial correlation of these antenna pairs directly to the 
source brightness distribution by a Fourier transform (Thompson et ah, 2008). Hence for radio imaging, 
each antenna pair output is cross-correlated to measure the coherence function which is subsequently 
converted to a sky map, conventionally by an inverse Fourier transform. Figure 1(b) shows the normalized 
Point Spread Function (PSF) corresponding to the aperture coverage in Figure 1(a), for a single point 
source along the w direction. The more densely sampled the projected aperture plane, the lower the 
spatial side-lobes of the sky image. The filling factor of the synthesized aperture can be increased by 
either using bandwidth synthesis or by populating sufficient baselines. In bandwidth synthesis different 
frequency channels can be used to scale A. As shown in Figurel(c) and Figurel(d), using only 10 frequency 
bins uniformly distributed across 1 — lOMHz, the aperture filling and the PSF is significantly improved 
as compared to Figurel(b). A first-order simulation of an array of V = 9 satellites in Earth-leading 
orbit around the sun yields Figurel(e) and the corresponding PSF in Figurel(f), where one snapshot 
each day is assumed at a single observation frequency of lOMHz. The number of uvw points are directly 
related to the unique number of baselines and the observational frequency. To achieve the confusion 
limit and resolve the sources individually, the total number of unique uvw points over the observational 
time period must be larger than the total number of detected sources. 
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(a) uvw: Snapshot (c) uvw: Bandwidth synthesis (e) uvw: Complete orbit 





Ou [arcmin] Ou [arcmin] Ou [arcmin] 

Figure 1: Baseline and Sky map simulations: Aperture filling of a 9-satellite ULW array for an Earth leading orbit 
around the Sun, to illustrate the effect of the sampling space on the normalized Point Spread Function (PSF). The uvw 
coverage for the 3-D array of satellites at at = lOMHz (a) for a single snapshot Nt = I, (b) single snapshot using 10 
frequency bins uniformly distributed in the range 1-10 MHz and (c) for an entire orbit around the sun at lOMHz with 
single observation each day i.e., 365 snapshots. 




2.2 Ultra-long wavelength sky 

The dominant foreground in the low frequency radio sky is the galactic synchrotron radiation, which is 
caused due to synchrotron emission from electrons moving in the Galactic magnetic field. This emission 
causes the brightness temperature to rise from ~10^K at 30MHz, to as high as ~10^K around 2MHz 
(Oberoi and Pincon, 2005). At frequencies below 2MHz, the Galactic plane is nearly completely opaque 
and the extra-galactic sources cannot be observed. More explicitly, for frequencies above 2MHz, the sky 
temperature can be approximated as (Jester and Falcke, 2009) 


r,.. = 16.3xI0«K(^) 


-2.53 


at > 2MHz, 


(3) 


where u is the observation frequency. For Earth-based observations at higher frequencies (> lOOMHz), 
the overall system noise temperature Tgys plaguing the cosmic signal is typically dominated by the 
noise from receiver electronics Tree- However, at lower frequencies (< 30 MHz), the intense galactic 
background implies that Tsky will be at least an order magnitude larger than T^ec, and hence the overall 
noise temperature Tsky ^ Tsys- The immediate effect of this extremely high sky noise is the poor 
sensitivity of the interferometric array. The 1-a RMS sensitivity for an antenna array of N nodes is 
(Cohen, 2004) 




_235.6 Tsys _ 

X^^/N{N - l)(uyihour)(Az//lMHz) 


mJy/beam, 


(4) 


where Au is the bandwidth, tobs is the observation time period over which the signal is integrated and 
the total number of estimated sources above this sensitivity is given by 


Ar>(S') = 1800(S'/10inJy) (i//10MHz) 


(5) 


Furthermore, the scattering in the interplanetary media (IPM) and interstellar media (ISM) also hinder 
observational frequencies less than 30MHz, which limit the maximum baseline between the satellites to 


disM = 47km x (z//lMHz)^'^ and diPM ~ 10km x (z^/lMHz). (6) 


Finally, the lower limit of the achievable noise is not the RMS sensitivity of the array, but the confusion 
limit. The presence of unresolved sources with individual flux densities below the detection limit leads 
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Parameter 

Notation 

Units 

Equation 

1 Extra-galactic survey | 

Sensitivity 


Jy 

Input 

6.5E-02 

6.5E-02 

6.5E-02 

6.5E-02 

6.5E-02 

6.5E-02 

Baseline 

d 

km 

Input 

100 

100 

100 

100 

100 

100 

Obs. Time 

^obs 

hours 

Input 

24 

720 

8760 

720 

8760 

8760 

Obs. frequency 


MHz 

Input 

10 

10 

10 

10 

10 

30 

Bandwidth 

diy 

MHz 

Input 

1 

1 

1 

3 

3 

3 

Resolution 

9 

arcmin 

(2) 

1.03 

1.03 

1.03 

1.03 

1.03 

0.34 

System temperature^ 

Tsys 

K 

(3) 

2.8E+05 

2.8E+05 

2.8E+05 

2.8E+05 

2.8E+05 

1.7E+04 

No. of Antennas 

N 


(4) 

229 

42 

12 

25 

7 

4 

ISM Max. Baseline 

diSM 

km 

(6) 

7.4E+02 

7.4E+05 

7.4E+05 

7.4E+05 

7.4E+05 

2.8E+06 

IPM Max. Baseline 

diPM 

km 

(6) 

100.0 

100.0 

100.0 

100.0 

100.0 

300.0 

Confusion limit 

^conf 

Jy 

(7) 

0.07 

0.07 

0.07 

0.07 

0.07 

0.01 

Resolution (Conf. lim.) 


arcmin 

(4), (7) 

1.00 

1.00 

1.00 

1.00 

1.00 

1.65 

Max. Baseline (Conf. lim.) 


km 

(2), (7) 

103.09 

103.09 

103.09 

103.09 

103.09 

20.86 

Time to Conf. lim. 


hours 

(8), (7) 

0.14 

4.05 

46.39 

3.98 

43.64 

38.94 


Table 2: System requirements for Extra-galactic surveys: The desired system parameters to achieve 
the desired resolution of 1' and sensitivity of 65mJy for Extra-galactic surveys, for varying observation time, observation 
frequency and instantaneous bandwidth. ( ^ Sky noise dominated) 


to a constant noise floor, that is reached after a certain observation time tobs (Jester and Falcke, 2009). 
For extra-galactic observations, under certain nominal assumptions, this anticipated confusion limit due 
to background sources is 

SconfiO, V) = 16mJy x (0eon//lT-"^(*^/74MHz)-O-7, (7) 

where 0 is the effective resolution for which the flux is below the confusion limit. The confusion limit 
is the lower limit to the achievable noise floor and thus is an upper limit to the useful collective area of 
the array. In other words, adding more antennas only decreases the time in which the confusion limit is 
reached, but not the overall array sensitivity (4). The time necessary for an array to reach achieve this 
confusion limited sensitivity is given by the “survey equation” 

tsurvey = 3.3 (Ar/lOO)-^ (lOi^/Aiv) (i^/lMHz) (8) 

Using these elementary and yet fundamental equations, a preliminary design for a space-based array can 
be proposed for desired science cases. For a more detailed study, refer to (Jester and Falcke, 2009). 

2.3 System definition 

The science cases for an ULW array broadly span cosmology, galactic surveys, transients from solar and 
planetary bursts and even the study of Ultra-High Energy particles. Although a single satellite mission 
would suffice to detect the global dark-ages signal, over 10^ antennas are required to investigate the radio 
emission from Extrasolar planets (see Jester and Ealcke, 2009, Table 1). Eor a first space-based ULW 
array however, with possibly only a few satellite nodes, extra-galactic surveys and study of transients 
are among the best suited science cases (Boonstra et ah, 2011), which we present as case studies. The 
proposed space-based array design can be readily extended to cater to other science cases e.g., detection 
of global dark-ages signal. 

The expected signal strength for the extra-galactic surveys are in the order of 65mJy with a desired 
spatial resolution of ~1'. In Table 2, we present different scenarios to investigate the effects of varying 
observational frequencies, bandwidth and observation time, on the number of antennas to achieve 65mJy. 
It is evident that increasing the observation time (1 day, 1 month, 1 year) steadily reduces the required 
number of antennas. Secondly, the increasing the bandwidth (IMHz to 3 MHz) is also an alternative to 
achieve the desired resolution for a small array. However, the increase in bandwidth has little effect on the 
confusion limit. We note that the confusion limit is a bottleneck for shorter integration times and lower 
observing frequencies. The maximum baseline is in general confusion limited for frequencies > lOMHz, 
however at < lOMHz, the ISM and IPM scattering limits the maximum baseline and subsequently the 
resolution. At lOMHz, we require at least one year of observation time with more than 7 antennas for an 
observational bandwidth of 3MHz to achieve the 65mJy sensitivity. However, in the last column of the 
Table 2, we see that at 30MHz, only 4 antennas sufficient. Such a configuration is estimated to detect 
over a million sources using (5). 

A similar investigation was conducted for Jupiter-like flares and Giant crab-like pulses, which are 
bright events with order of MJy and kJy respectively, with typical time scales of milliseconds. The 
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Figure 2: Space-based antenna model, configurations and simulations: (a) System model for an 
LNA connected to the Antenna (b) Configuration of a dipole antenna on a cubesat (c) Configuration of a tripole antenna 
in a cubesat (d) four monopoles (e) A dipole antenna placed on the opposite ends of a 3-cubesat (f) Required as a 

function of frequency for non-matched dipoles of varying lengths, (g) Real part of input impedance and (h) Absolute value 
of imaginary part of input impedance for two monopoles placed on a satellite compared to a dipole (i) Ratio of gain of two 
monopoles placed on a satellite and a dipole. 


desired resolution for these transient radio systems are at < 30MHz. Since these events are extremely 
bright, even a single antenna with a nominal bandwidth of 10 % of observational frequency would meet 
the desired requirements. However, unlike the extragalactic surveys the observations are not confusion 
limited but possibly by the number of baselines for short integration timescales of milliseconds. The 
number of unique uvw points will depend inadvertently on the deployment location and the relative 
range rate of the antennas. However, in general this limitation can be overcome by increasing the 
integration times in both these cases by over a minute. 

In general, higher bandwidth, higher observing frequencies and longer integration times require less 
antennas to reach the same sensitivity level. In this article we choose the DARIS mission specifications 
as a reference to illustrate various sub-systems. To this end, we particularly focus on an array of 
N = 9 satellites, with a maximum satellite separation of 100 km and capable of observing the skies at 
O.I — lOMHz. This particular setup meets the requirements for the extra-galactic survey and transient 
radio system science cases. However, all the proposed techniques and technologies can be readily extended 
to a larger array, for observation frequencies up to 30MHz. 

3 Radio astronomy antenna design 

A critical component for the space-based array is the design of the observational antenna. The system 
model of the observation antenna connected to a LNA is shown in Figure 2(a), wherein the antenna is 
modeled as an ideal lossless antenna, followed by an attenuator representative of the antenna losses. For 
the observation frequencies of 0.3 — 30MHz, this front end must be sky noise limited i.e., T^ec < 
where Tree is the receiver noise and Tsky is the sky noise temperatures which are defined at the interface 
between the lossless antenna and attenuator. The LNA noise temperature Tlna defined at the input of 
the LNA is equal to (I — 7 ^)To, where 77 is the radiation efficiency and Tq is the physical temperature of the 
antenna (chosen as 290K). Without loss of generality, we assume that the LNA noise is dominant over the 
noise contribution of subsequent electronics of the receiver. Under these assumptions, the prerequisite 
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on the LNA noise temperature is derived as 


Tlna < (1 - |r|^)(0.l7^Ts/c^ “ (1 ~ (9) 

where F is the reflection coefficient of the antenna (Arts et ah, 2010). A straightforward candidate 
for the observation antenna is a dipole (e.g., Figure2(b)), which can be realized by rolling out metallic 
strips from the satellite (Manning, 2000). The observational wavelengths are much larger compared to 
the dimensions of the satellites and hence due to practical limitations, the realized dipole will be short 
compared to the wavelength. For instance, a classic half-wave dipole for lOMHz and 30MHz observation 
frequencies yields a dipole length of 15m and 10m respectively. For lower frequencies, a similar dipole 
lengths begets a short-dipole with a directional pattern similar but less directional as compared to the 
half-wave dipole. Consequentially, the radiation resistance of the small antenna would be low and the 
thermal noise will significantly dominate the total antenna noise. 

To maximize the received power at the antenna, matched dipoles can be used, in which case F ^ 0. 
However, the combination of antenna and matching network becomes highly resonant with a high quality 
factor, significantly limiting the achievable bandwidth of the system, in particular for shorter dipoles. 
Hence, we propose the use of a non-matched dipole (Arts et ah, 2010). Using (9), we have the Figure2(f), 
which shows the Tlna for varying lengths of non-matched dipole lengths. We use a LNA with input 
impedance of 2KU, which is a realistic value for an LNA using a bipolar transistor. While the Tlna 
is above lOOK for 10m and 15m antennas, the required noise temperature for shorter dipole lengths are 
significantly lower, especially for lower frequencies. 

In practice, a dipole is implemented on a satellite using two monopoles. To verify the validity of the 
proposed model, impedances of two monopoles on a satellite body is compared to that of the dipole. 
The satellite body under simulation is modeled as a cube of 40 x 40 x 40cm , with perfect conducting 
surfaces. Furthermore, since the impedance of the monopole above a perfect ground plane is half the 
impedance of a dipole, the impedance of the monopole is multiplied by a factor 2 in the simulations for 
a fair comparison. Figure 2(g) and Figure 2(h) show the real component of the impedance and absolute 
value of the imaginary component of the impedance respectively, for Im dipole versus 0.5m monopoles 
and 15m dipole versus 7.5m monopoles respectively. As seen in these figures, there is negligible difference 
between the dipole and the two monopole configuration. A step further, we compare the ratio of gains 
between a pair of monopoles against the 1 meter and 15 meter dipoles at 1 and lOMHz. The investigated 
antenna lengths the ratio of gain are almost flat across for varying angles as seen in Figure 2(i), which 
indicates the element pattern does not change if a configuration of two monopoles is used instead of a 
dipole. 

Two orthogonal dipoles (Figure2(b)) are in principle sufficient to get all the polarization information 
of the cosmic signal, however a tripole (i.e., three dipoles, see Figure2(c)) can be used to obtain infor¬ 
mation of all 3 components. The third dipole improves the directivity of the antenna system, thereby 
increasing the field of view. Along similar direction, a equiangular spaced four monopole configuration 
can also be considered, as shown in Figure2(d). However, the number of correlations is much higher 
and consequentially demanding more signal processing hardware for each antenna. In the OLFAR study 
where a 31/-cubesat (30 x 10 x 10 cm) is utilized, the monopoles are deployed in groups of three at the 
opposite ends of the satellite, as seen in Figure2(e). The asymmetric design changes the properties of 
the monopoles and reduces the purity of the independent components, which is studied by Smith and 
Arts (2013) and experimentally evaluated on a smaller scale by Quillien et al. (2013). 

4 Digital Signal Processing 

Sky images in radio astronomy are made by calculating the fourier transform of the measured coherence 
function (Taylor et ah, 1999). The coherence function is the cross correlation product between two 
antenna signals located at the two spatial positions, averaged over a period of the integration time Tint- 
One way to implement such a system is using the traditional XF correlator Le., cross correlation first 
and fourier transform later and the more recent FX correlator which measures directly the cross-power 
spectrum between the two antenna signals (Bunton, 2004). Although the XF architecture is beneficial 
because bandwidth can be traded for spectral resolution, FX architecture offers computational efficiency 





Figure 3: Correlator architectures: An illustration of two potential correlator architectures for space-based radio 
interferometric array, where the tags ‘X’ and ‘T’ on the nodes indicate Correlation and Transmission to Earth capabilities 
respectively. In the (a) Centralized correlator architecture a centralized mother ship receives data from all observational 
satellites, correlates and down-links data to Earth. On the contrary, in the (b) Distributed correlator framework, the 
observed data is evenly distributed between all nodes. Post correlation, all satellite nodes down-link their respective 
correlated data to Earth. 


and more significantly. The processing factor for XF vs FX is given by 

J^xflfx _ / NsigNbins \ 

^ ^sig T ^hins ' 


( 10 ) 


where Ngig = NpoiN and N^ins = ^fi/^fres (Rajan et ah, 2013b). Observe that the multiplicands 
in the XF mode are additive in the FX mode , besides the log2 reduction on the number of frequency 
bins. Thus, although for lower number of nodes the XF is comparable to FX mode, for large scalable 
architectures the FX mode is computationally cost effective. Since we prefer a scalable space based array 
the FX architecture is chosen as the preferred architecture . Table 3 shows data rates for cluster of 
N = 9 nodes, with instantaneous bandwidth of Afi = IMHz and r = 1 second integration time. 

A typical pre-processing unit at each satellite node is shown in Figure 4, where each satellite gen¬ 
erates Dobs = ‘^AfiNpoiNbits^ps. Observe that with Npoi = 3, for a signal with nominal instantaneous 
bandwidth of Afi = IMHz sampled with Nbu = Ibit resolution, the output data rate is 6Mbps. Given a 
far away deployment location, such as Lunar orbit (~400, 000 km) or Earth leading/trailing (~2 x 10^ to 
~4 X 10^ km), this down-link data rate levies heavy prerequisites on the limited resources of a small satel¬ 
lite using current techonology. Hence, the satellite cluster must not only employ on-board pre-processing 
of astronomical signals, but also on-board correlation to minimize down-link data rate back to Earth. 
To this end, either a centralized or a distributed correlator can be employed as illustrated in Eigure3. 


4.1 Centralized architecture 

In the centralized EX correlator framework each satellite node transmits Dobs = ‘^^fi^poi^bits^P^ 
to the centralized mothership, which in turn receives D^^ = Dobsi^ — l)bps in total, excluding the 
data collected from the antenna on the mothership itself. The input data from all satellites is then 
correlated and the output is then transmitted down to Earth at the rate D^^^ = 2N‘^^gNbins/ 
where Nbins = ^fi/^fres- A significant drawback of the centralized correlation is that it depends 
heavily on the healthy operation of a single correlation station, which introduces a Single Point Of 
Eailure (SPOE) for large array of satellites in space. 

4.2 Distributed architecture 

To alleviate SPOE a Frequency distributed correlator is proposed where each node is pre-assigned a 
specific sub-band Afsb of the observed instantaneous bandwidth Afsb for cross correlations (Rajan et ah, 
2013b). Hence, in addition to the node pre-processing (Eigure4), a secondary fine PEB is implemented to 
further split the instantaneous bandwidth Afi into Nsb sub bands, each of bandwidth Afsb = Afi/Nsb- 
Each satellite is assigned a specific sub-band for processing and the other {Nsb — 1) sub-bands are 
transmitted to corresponding satellites via the intra-satellite communication layer. Eurthermore, for 
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Figure 4: Node level signal processing: Given the low observational bandwidth of A/o < 50 MHz, the Npol = 3 
polarized astronomical signals received by each antenna will be signal conditioned and Nyquist sampled at 2A/o MHz with 
a 14-bit (or more) Analog to Digital Converter (ADC) as shown in Figured. A coarse Poly-phase Filter Bank (PFB) is 
used to selectively choose the desired instantaneous bandwidth of Afi = IMHz. After successful RFI Mitigation (RFIM), 
only = 1 — 2bits will be used in further processing stages (Altunin, 2001). The total data generated for Np^i signal 

paths in each satellite is = OA/^A^^^gbps, which is transported to the Inter-satellite communication layer. 


Data rates and processing 

Notation/Equation 

Value 

U nits/Remark 

No. of satellites (or antennas) 
No. of polarizations 

No. of channels/signals 

No. of bits 

Observation frequency range 
Instantaneous bandwidth 
Spectral resolution 

No. of bins 

Snapshot integration time 
Observed data rate 

N 

Ftpoi 

Ff sig — FfpolN 

F'bits 

Ffo 

A/. 

Ffres 

Fbins — Ffij XfpQs 

T 

Cobs — ‘^FfiNpolNbits 

9 

3 

27 

1 

-^30 

1 

1 

1000 

1 

6 

(scalable) 

bits 

MHz 

MHz 

kHz 

second 

Mbps/satellite 

Centralized 




mothership data reception 
Earth down link data rate 

= DoUN - 1) 

Cout ~ ‘^Fg^^NbitsFbins 

48 

1.46 

Mbps 

Mbps 

Distributed 




No. of sub-bands 

Sub-band bandwidth 
Inter-satellite reception 

Earth down link data rate 

Nsb=N 

Afsb = Afi/Nsb 

Dt = DUN 

Dint = DU/N 

9 

111.11 

5.34 

162.2 

kHz 

Mbps/satellite 

kbps/satellite 


Table 3: Centralized vs Distributed processing: Data rate estimates for a Centralized eorrelator 
and a Frequeney distributed FX eorrelator for the DARIS mission of 9-satellites. 


even distribution of data and to ensure scalability, we enforce the number of sub-bands equal to the 
number of satellite nodes Ae., Nsh = N. Subsequently, in the distributed framework, each node receives 
a specific sub-band of the observed data i.e., {Dobs/^sh) from N — 1 other satellites in the network which 
yields a total input of = {Dobs/Nsb){N — 1) = (D^^/A^)bps, and down-links = {D^p^^/N)hps 
respectively. 

Thus, the Frequeney distributed eorrelator reduces the inter-satellite communication by a factor N. 
Furthermore, at the cost of equipped quipping all observational satellite nodes with communication 
capability (both Inter-satellite and down-link to Earth), SPOF is averted and scalability is ensured. In 
the context of the projects discussed earlier, DARIS, FIRST and SURO-LC implement a centralized 
architecture, whereas OLFAR employed the distributed architecture. Given that the system frequency 
is typically an order or more than the processing instantaneous bandwidth, computing requirements are 
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Short term clock stability 



Figure 5: Clock stability : (Left) Short term: The plot shows limiting cases of the Signal to Noise Ratio (SNR) and 
corresponding Effective Number of Bits (ENOB) due to jitter tjitter versus input frequency Uin . Three demarkation lines 
shows the maximum input frequency of Uin = lOMHz, Uin = 30MHz and Uin = 50MHz. (Right) Long term: Desired Allan 
deviations of free running clocks are plotted versus the coherence time (in green) for various input frequencies Uin- The 
map is overlayed with Allan deviations of potential clocks (in blue) for potential space-based low frequency arrays namely 
PRS-10 Rubidium (sta, 2006) , RAFS ASTRIUM (Droz et ah, 2007), GPS Ipps (Lombardi et ah, 2001), OCXO ASTRIUM 
(ocx, 2012), SA.45s CSAC and Space Hydrogen Maser (SHM) (Goujon et ah, 2010). 


negligibly small, which has been duly noted by all these studies. 

4.3 Clocks 

The choice of the on-board clock on each satellite has a significant impact on the signal processing 
system . The short-term clock stability i.e., t 1 second is dominated by the clock jitter, which limits 
the Effective Number Of Bits (ENOB) for a given input frequency Uin. Eor instance, to facilitate 14- 
bit sampling at Uin the chosen clock must have dtjuter < Ips, as seen in the Eigure 5(a). Secondly, 
to evaluate clock stability over long durations Ae., t ^ 1 seconds, we use Allan variance^ which is a 
measure of nominal fractional frequency drift (Barnes et ah, 1971). Eollowing (Thompson et ah, 1994; 
Ulvestad et ah, 1986), the stability requirement on the clock and define the coherence time Tc, such that 
the RMS phase error of the clock remains less than 1 radian 

iy(T(^{Tc)rc ^ 1 , ( 11 ) 

where n is the observational frequency and (J(^{tc) is the Allan deviation as a function of Tc (Rajan 
et ah, 2013a). The product (T(^{tc)Tc can be visualized as the time drift due to non-linear components 
of the clock after Tc seconds. Eurthermore, the linear parameters of the clock i.e., frequency and phase 
offsets can be eliminated by exploiting the affine clock model, which is discussed in Section 6. Eigure5 
shows expected Allan deviations of potential clocks versus the coherence time as per (11) for various 
input frequencies vin. Among the presented choices in Eigure5, the RAES ASTRIUM (3.3kg, 30W) and 
OCXO ASTRIUM (220g, 2W) are space qualified Rubidium and Oven controlled oscillators respectively. 
A particular clock of interest is the SA.45s, which is a Rubidium clock weighing less than 35 grams, 
consuming < 0.125 watts and offers an coherence time of up to 15 minutes. 

5 Communications 

The potential communication scenarios for the envisioned space-based array are shown in EigureG, which 
follow directly from the correlator architectures discussed in the previous section. A centralized architec¬ 
ture, as shown in EigureG (a), comprises of a mothership collecting raw observed data from a cluster of 
daughter satellites and, down-links the processed data to an Earth-based ground station. Alternatively, 
in case of the distributed scenario shown in EigureG (b), all satellites are capable of both exchanging data 
and correlating them, before down-linking back to Earth. In addition to the science data, housekeeping 
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Figure 6: CoTTHnUTlicCitioTl architectures: An illustration of a (a) Centralized communication architecture 
and (b) single pairwise-link of a Distributed communication architecture for a space-based radio interferometric array. 
The Inter-satellite link is indicated in blue and the Earth-downlink by red . Telemetry and Tele-commands are exchanged 
between the satellites and with Earth in both scenarios. 


information is also exchanged between the satellites via tele-commands and telemetry, which is expected 
to be relatively small 100kbps) in comparison to the astronomical data of 6Mbps. The housekeeping 
information is critical for control, timing and synchronization of the satellite, and to maintain coherence 
within the satellite network. 

5.1 Inter-Satellite Link (ISL) 

Implementing the ISL using high-frequency optical communication (Sodnik et ah, 2006; Toyoshima, 2005) 
has many advantages as compared to radio communication, such as reduced mass and volume of equip¬ 
ment, higher data rates and no regulatory restrictions as experienced for Radio Frequency (RF) bands. 
However, this would also require extremely accurate alignment of the satellites, robust synchronization 
and more power than what could potentially be available for a small satellite. In the RF domain, OFDM 
is an efficient modulation scheme for the ISL, in particular for a scalable antenna array with limited 
available bandwidth (Nee and Prasad, 2000). The main advantages of OFDM are it is very well suited 
to frequency selective channels and it potentially offers a good spectral efficiency. With OFDM, the sep¬ 
aration between each channel is equal to the bandwidth of each channel, which is the minimum distance 
by which the channels can be separated. The signals from each satellite node which form individual 
channels will be modulated using a form of Phased Shift Keying-PSK, Amplitude Shift Keying-ASK, or 
a combination Quadrature Amplitude Modulation-QAM. In this article, we consider an ISL transmission 
frequency of 2.45GHz, although other frequency bands can also be used. 

One of the possible solutions to implement the ISL is to use patch antennas on each face of a satellite 
node, such that the combined implementation yields a full coverage of the sky. All the satellites will 
have patch antennas on all six faces for both uplink and downlink. In addition, a diplexer will be used 
to separate the receiving and transmitting channels. Using a coaxial switch (controlled by the CPU) 
this signal is selected, amplified and finally detected in the node. The desired antenna must have a 
bandwidth of lOOMHz around 2.4GHz with the reflection coefficient less than —10 dB between 2.35GHz 
and 2.45GHz, and the half-power beam-width must be at least 90°. Figure 7(a) shows the simulated 
patch antenna, where all dimensions are in millimeters. The patch is fed by a coaxial probe located 
7.5mm from the center of the patch. It was observed that in case of a linearly polarised patch antenna, 
the given setup yields a —3 dB beamwidth less than 90° for the radiation pattern in the (j) = 90° plane, 
which is less than our desired requirement. Thus, to improve the radiation pattern we implemented 
a circularly polarized patch antenna. Moreover, an added advantage is that the polarization of both 
transmit and receive antenna is independent of the orientation of the antennas with respect to each 
other. The circular polarization is realized by adding a second co-axial probe to the patch, shown in 
red in Figure 7(a). As seen in the Figure 7(b), the reflection coefficient is better than —10 dB in the 
frequency range 2.35 — 2.8 GHz. Additionally, the radiation patterns in Figure 7(c) show the half-power 
beam width is more than 3dB for both the 0 = 0° and 0 = 90° planes. 
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Figure 7: Patch AntCTlTia for Inter-Satellite Link (ISL): (a) The simulated patch antenna with dimensions 
in millimeters, (b) Reflection coefficient and (c) Radiation patterns across desired frequency range. 


5.2 ISL Link margin 

Table 4 shows the ISL Link budget for the centralized and distributed scenarios. In case of the centralized 
scenario, we assume that the mothership is in the center of the array with a maximum mothership-Node 
distance of 50km. To estimate the link margin for a centralized scenario, we assume that the mothership 
can transmit with a power of IW i.e., 30dBm. Now, using the patch antenna gain of 3dBi and nominal 
losses (at the diplexer, transmitting cable and coaxial switch loss) of 2.6 dB, the EIRP (Equivalent 
Isotropically Radiated Power) of each satellite is 30.4dBm. Eor both centralized and distributed scenarios, 
the ISL channel in space is in principle free space loss, where Multi-path, atmospheric losses, absorption 
losses and even Doppler effects can be ignored. Hence, the Eriis free space loss for 2.4GHz transmission 
frequency and 50km is 134.7dB. Hence, the effective received power for an EIRP of 30.4dB and a pointing 
loss of 0.5dB is —104.80dB. The received C/NO is estimated at 67.66dB/Hz, for a G/T ratio of —26.14 
dB/K (see Table 4). The transmission data rate mothership to Node is 100kbps (50 dB/Hz), which 
yields an Eb/No of 17.66dB. Eor a typical receiver, an Eb/No of 2.5dB is needed. Now, including a 
implementation loss of 2dB, the link margin for uplink with transmit power of IW is 13.16dB. The 
downlink from the Node satellite to the mothership, includes the 6Mbps data (see Table 3) and the 
housekeeping data of 100kbps, which amounts to 6.10Mbps. This downlink can be achieved with a link 
margin of 2.29 with a transmission power of 5W. 

Extending the link margin estimates of the centralized ISL architecture to a distributed scenario 
has two fundamental challenges. Eirstly, the transmission data rate between every satellite pair is now 
5.44Mbps, which includes 5.34Mbps of science data (see Table 3) and 100Kbps of housekeeping data. 
Secondly, in the absence of a centralized correlator, the maximum distance between the satellites is 
100km, a factor 2 compared to the centralized scenario. Hence, to achieve the same link margin of 
2.29dB as the Node to mothership downlink, the transmission power of each satellite in the distributed 
architecture must be 4 times that of the centralized scenario, i.e., 20W. Although, 15W suffices to achieve 
a positive link margin for the distributed architecture. This requirement is a bottle neck for scalable 
array of small satellites with limited transmission power. One possibility is to use clustering schemes 
and multi-hop approaches to reduces the communication distances between the satellites (Budianu et ah, 
2011), which is a research area currently being explored (Naghshvar et ah, 2012). 

5.3 Space to Earth Downlink 

The total downlink data rate Dout after correlation quadratically increases with the number of nodes in 
the cluster and reduces linearly with the integration time (see Table 3). Eor a cluster of 9 satellites, with 
1 second integration interval this rate is > 1.46Mbps. In the DARIS study, the space to earth downlink 
was achieved using an X-band Downlink Unit (XDU). Eigure 8 shows the X-Band transmit chain which 
consists of a modulator and a 60W Traveling Wave Tube Amplifiers (TWTA) with 60% efficiency. In 
the nominal operation case the data is directly BPSK modulated on the X-band carrier at 4.3Mbps, 
amplified, filtered at the diplexer and transmitted via the high gain antenna (40dBi). During launch and 
early orbit phase and in contingency cases the tranmission rate is reduced to 100bps, BPSK-modulated 
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Parameters 

Centralized 

MS Node Node MS 

Distributed 

Node Node 

Units 

System Definition 

Frequency band 

2.45 

2.45 

2.45 

GHz 


Maximum distance 

50.00 

50.00 

100.00 

km 


TX power 

1.00 

5.00 

15.00 

W 


Required datarate 

0.10 

6.10 

5.44 

Mbps 


Antenna gain 

3.00 

3.00 

3.00 

dBi 


Wavelength 

0.12 

0.12 

0.12 

m 

ISL Channel 

Free space Path losses 

134.20 

134.20 

140.23 

dB 


Atmospheric losses 

0.00 

0.00 

0.00 

dB 


Polarization losses 

0.50 

0.50 

0.50 

dB 


Absorption losses 

0.00 

0.00 

0.00 

dB 


Path loss 

134.70 

134.70 

140.73 

dB 

Link Budget 

TX power 

30.00 

36.99 

41.76 

dBm 


Loss coaxial switch 

0.30 

0.30 

0.30 

dB 


Diplexer losses 

1.00 

1.00 

1.00 

dB 


Cable losses 

1.30 

1.30 

1.30 

dB 


Transmit antenna gain 

3.00 

3.00 

3.00 

dBi 


EIRP of the spacecraft 

30.40 

37.39 

42.16 

dBm 


Pointing loss 

0.5 

0.5 

0.5 

dB 


Received power 

-104.80 

-97.81 

-99.06 

dBm 


Receiver G/T 

-26.14 

-26.14 

-26.14 

dB/K 


Bolzmann constant 

-198.60 

-198.60 

-198.60 

dBm/Hz/K 


Received C/NO 

67.66 

74.65 

73.40 

dB/Hz 


Data rate 

50.00 

67.85 

67.36 

dB/Hz 


Received Eb/No 

17.66 

6.79 

6.04 

dB 


Required Eb/No 

2.5 

2.5 

2.5 

dB 


Implementation loss 

2 

2 

2 

dB 


Link Margin 

13.16 

2.29 

1.54 

dB 


G/T estimates 

Gain (dB) 

Temp (K) 

Temp (dBK) 

Loss (dB) 

Tref (dBK) 

Gref (dB) 

Tref (K) 

Antenna 

3.00 

100.00 

20.00 

-2.40 

17.60 

3.00 

57.54 

Coaxial Cable 

-1.30 

75.02 

18.75 

-1.10 

17.65 

-1.30 

58.23 

Diplexer 

-1.00 

59.64 

17.76 

-0.10 

17.66 

-1.00 

58.29 

Coaxial Cable 

-0.10 

6.60 

8.20 

0.00 

8.20 

-0.10 

6.60 

LNA 

30.00 

288.63 

24.60 

0.00 

24.60 


288.63 

Backend 

30.00 

2610.00 

34.17 

-30.00 

4.17 


2.61 

LNA Noise Figure (dB) 

3.00 




G (dB) 

0.60 


Backend Noise Figure (dB) 

10.00 



Totals 

T (K) 


471.90 

Ambient Temperature (K) 

290.00 




G/T (dB/K) 


-26.13 


Table 4: Inter-sallite communication: (Above) The Inter-Satellite Link (ISL) budget for centralized 
and distributed scenarios for a satellite array. (Below) Antenna gain to noise tempearature (G/T) 
estimates. 


onto a subcarrier with low frequency (e.g., 8kHz) which itself PM-modulates the carrier. In order to 
enable full coverage two adversely polarized low gain antennas are transmitting into the two hemispheres. 
The receiver downconverts and demodulates the X-band data and feeds them to the OBC via RS422 link. 
In the nominal mode the command rate is 400kbps, BPSK-modulated on a 16kHz subcarrier and PM 
modulated onto the carrier. The estimated power consumption of the transceiver is IlOW operational 
and SOW in standby. The HGA, covering both up- and downlink band, could be realized as a parabolic 
antenna similar to the one of Mars Express. With a 40dBi design, the 3dB beamwidth would be around 
1.7° posing no big challenge to the pointing mechanism. The LG A is a conventional helix antenna design 
covering also both transmit and receive band. The number of ground stations on Earth are limited 
and are therefore almost always in high demand. Eor instance, the core ESA network has two deep 
space network 35m antennas along with several antennas in the 13 — 15m class (Vassallo et ah, 2007). 
A baseline for ground station access for is ESAs 35m ground station for 8hours/day suffices the need 
for the data generated cluster of 10 satellites with a centralized mothership. Eor arrays larger than 50 
satellites with minimal power, it is difficult to establish Earth-based downlink with current off the shelf 
technology. However, these challenges and possible distributed downlink scenarios are currently being 
investigated (Budianu et ah, 2015). 

6 Synchronization and Attitude determination 

To maintain coherence between satellites, all the satellites must be synchronized, and their positions 
known up to sub-meter accuracies. It is worth noting that these requirements on space-time accuracies 
are much lower compared to other space-based array missions, e.g., LISA (Bik et ah, 2007). Almost all 
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Figure 8: DARIS mother ship-Earth Link X-band Downlink Unit (XDU): The XDU unit consists 

of a complete redundant transmit and receive subsystem, the X-band antennas and the interconnecting waveguide links. 
The X-Band downlink transceivers are directly interfaced to the Onboard Computer (OBC) via RS422. Two different 
modes are available for the communication link (a)Nominal link via High Gain Antenna (HGA) for science telemetry and 
commanding (b) Emergency link via low gain antennas (LGA) with full coverage and reduced data rates for telemetry and 
command. 


Earth-based antenna arrays synchronize using GPS-aided atomic clocks, where fixed antenna positions 
are known up to millimeter accuracies, cf. e.g., LOFAR (van Haarlem et ah, 2013). However, the 
envisioned space-based array will be deployed far-away from Earth-based GPS satellites and unlike Earth- 
based antennas, these satellites will be mobile within a formation or swarm due to orbital dynamics. 
In addition, given the large number of satellites and limited ground-segment capability, tracking each 
satellite independently is infeasible. Moreover, in certain deployment locations such as the lunar orbit, 
the satellite array will be partially or even completely disconnected from Earth-based ground stations 
during eclipse. 

Hence, the satellite array must be a cooperative wireless network^ collaboratively estimating time- 
varying positions and correct for respective clocks simultaneously. More specifically, the satellites will 
estimate relative positions which suffice for Radio astronomy imaging, inter-satellite communication and 
collision avoidance. The relative positions can be estimated from distances via MDS-like algorithms 
(Borg and Groenen, 2005) (Rajan et ah, 2015). In this section we estimate the clock parameters and the 
time-varying distances of the network. We are primarily interested in solutions for cold start scenarios, 
when no prior information is known. For longer time scales, when the orbital dynamics of the deployment 
location is well known, the estimated space-time parameters can be tracked and improved using recursive 
filters e.g., Kalman Filter (Kay, 1993). The pointing direction for the satellites can be provided by 
commercially available sun (or star) trackers, which form part of the Attitude and Orbit Gontrol System 
(AOGS) in the satellites. 

6.1 Joint ranging and synchronization 

All clocks are inherently non-linear w.r.t. the ideal time t. However, a given clock can be approximated 
to a linear model provided the Allan-deviation of the clock is relatively low for a small coherence time 
(see Section 4.3). More generally, let U^tj denote the local times at i,jf respectively, then the ideal time 
t is 


t= C{ti, ipi) = C{tj, -Ipj, ipj), (12) 

where and are the phase and frequency offsets of a satellite pair (i^j) and, C(') represents 

a linear function of the local clock parameters. However, in case a satellite network, the pairwise distances 
are time-varying, which can be expanded as 


^ij (^) ~ '^ij ij^ ^ij^‘^ + • • • , (13) 

where dij (t) is the time-varying distance between the satellite pair (i^j) and {rij , gj, icij ,...) are range pa¬ 
rameters of the Taylor expansion at ideal time t = 0. More specifically, rij is the initial pairwise distance, 
rij is the range rate and Gj denotes the rate of range rate between the satellites. Given these pairwise 
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Figure 9: DyTlCLUlic VdTigiTig: A generalized Two Way Ranging (TWR) scenario between a pair of asynchronous 
mobile satellite nodes where the asynchronous nodes transmit and receive asymmetrically, during which K time stamps 
are recorded at respective nodes. Unlike classical TWR (iee, 2007) where the transmission and reception is alternating, 
the proposed setup imposes no pre-requisites on the sequence or number of two way communications. Consequently, this 
framework and proposed solutions for joint ranging and synchronization (Rajan and van der Veen, 2015) can be readily 
extended to a plethora of TWR ranging protocols, including broadcasting and passive listening 

distances, the relative positions of the satellites can be estimated using MDS-like algorithms (Borg and 
Groenen, 2005). Our aim is to jointly estimate the clock parameters and the time-varying pairwise dis¬ 
tances between the satellite nodes. The joint synchronization and ranging problem can be formulated 
as shown in Figure 9, which shows a pair of asynchronous mobile satellite nodes. The mobile satellites 
transmit messages asymmetrically between each other, during which K time-markers are recorded at 
each end. Let k and k be the kth time-markers recorded at the ith and jth satellite nodes respec¬ 
tively, and Eij^k C {+1,-1} indicates the transmit and receive direction of the message. For any kth 
time instant, the Generalized Two-Way Ranging (GTWR) equation (Rajan and van der Veen, 2015) is 

+ Eij^kdij(Ci{Tij^k,^i,'>Pi)^ = 0 Vz,j < N ,'ik < K, (14) 

where without loss of generality the ideal time t of the time-varying distance dij (t) is replaced with the 
clock model at satellite i (12). The unknown frequency offsets, phase offsets and the pairwise distances 
over K time instances can be estimated using the iterative Mobile Pairwise Least Squares (iMPLS) and 
iterative Mobile Global Least Squares(iMGLS). The iMPLS algorithm is applicable when the daughter 
satellite nodes communicate only with a centralized Mother-ship (i.e.. Star network topology, see Figure 3 
(a)), whereby only the clocks of the satellites can be corrected. For a full mesh network, when all satellites 
communicate with each other (i.e.. Full mesh topology, see Figure 3 (b)) both the clock parameters and 
the distances can be estimated, which is achieved by the iMGLS algorithm. 

6.2 Simulations 

To illustrate the algorithms, we consider a cluster of V = 9 mobile nodes in a 3—dimensional Euclidean 
space. The initial frequency offsets -0 = ['0i,'02, • • • ,'^ 9 ] and phase offsets -0 = ['0i,'02, • • • of the 
nodes are arbitrarily chosen in the range [—10“^, 10“^] and [—1,1] seconds respectively, where without 
loss of generality satellite node 1 is chosen as the reference node with = [0,0]. Secondly, the 

initial positions X and initial velocities Y, whose values are arbitrary chosen as 

" 6.1 -7.6 -0.2 8.2 9.0 -8.9 9.9 6.0 -9.3“ 

-1.3 4.4 2.9 -9.7 -1.8 -7.7 -7.6 7.8 7.8 km, 

_ 3.5 6.8 9.1 5.3 5.4 6.6 6.0 4.6 7.5 _ 

“-70 -50 -30 -70 -70 40 -40 60 70l 

90 40 40 -90 -100 10 -80 -40 70 ms“\ (15) 

_ 60 50 20 50 80 60 100 90 10_ 

For these assumed positions, the initial pairwise distance r is in the range [0,10]km. Secondly, for 
these values the range rates r are distributed in the period [—100,100]m/s and the rate of range rate 
r are distributed over [—10,10]m^/s. These values are typically for worst case scenarios, assumed to 
evaluate the performance of the proposed algorithms. The mobile satellites communicate K messages 


X = 

Y = 
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Frequency offset, '0 Phase offset, 0 



Figure 10: Joint ranging and synchronization: A simulation showing the Root Mean Square Errors 
(RMSEs) of the estimated frequency offset (0), phase offset (0), range parameters (r, r, r) and distance 
d using the MGLS algorithm for joint ranging and synchronization. 


with each other within a time interval of [—10, lOjseconds, and the time-markers K are 

generated accordingly. We assume alternating communication between the nodes and a Gaussian noise 
of cr = 10“^ seconds~3.3 meters plaguing the time-markers. Eigure 10 shows the performance of the 
proposed algorithms for phase offset, frequency offset and pairwise distance estimation. The Root Mean 
Square Errors are plotted (in blue) against varying iF = 10toiF=100, where the clock parameters 
are averaged over N nodes and, the distances and range parameters are averaged over (^) unique links. 
In addition, the Cramer Rao Lower Bounds for the corresponding estimates are also plotted (in red), 
which is achieved asymptotically by the proposed estimators. To show the performance of the prevalent 
solutions, we also plot the Low Complexity Least Squares (LCLS) which synchronizes the clocks for 
immobile network of satellite nodes. The iMGLS algorithm outperforms the iMPLS estimator since 
it exploits the full mesh communication network between the nodes. More significantly, the iMGLS 
estimator achieves clock accuracies upto nanoseconds and distance errors up to meter accuracies at cold 
start. 

An added advantage of using dynamic ranging is that the timestamps can potentially piggyback on 
the housekeeping data exchanged between the satellite nodes, which mitigates the need for a dedicated 
ranging system. However, if a ranging system is employed, then the achievable lower bound on the 
standard deviation for Time Of Arrival in multipath-free channels is given by 

O- > {Stt^F^BT SNR)“^^^ ( 16 ) 

where Fc denotes the carrier frequency, B Fc is the bandwidth of the signal, T is the signal duration 
in seconds (Patwari et ah, 2005) . The assumed noise variance on the time-markers in the simulation 
is cr = 10“^ (shown in green in Eigure 10), which can be adequately achieved by a wireless node 
communicating at Fc = 2.4GHz with a nominal bandwidth of IkHz transmitting and SNR=10dB for a 
signal duration of T~lms. 
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Figure 11: Orbiting the Moon: (Left) Relative orbits with different relative inclinations. (Right) Two relative 
orbits with several spacecraft in each orbit. These figures represent the relative orbits with respect to a reference orbit, 
with the Moon as the central body, without any perturbations. 


7 Deployment Locations 

The deployment location of the space-based array must be chosen to ensure the following conditions. 

• Minimize RFI during scientific observation cycles. 

• Offer maximum possible down-link data rate. 

• Provide sufficient positional stability during integration time r. 

• Must remain within a sphere of ~ 100km. 

In addition, each satellite must offer low noise conditions with minimal EMC and stable temperature (and 
gain) conditions to easen calibration easier. To alleviate the high complexity of active control to keep all 
the satellites within a cluster, passive formation flying could be employed. In passive formation flying 
paradigm, the satellites are allowed to drift, during which the relative positions and orientations of the 
satellites are constantly monitored. This approach eliminates the need for excess propulsion and heavy 
orbital maintenance equipment on all satellites. Additionally, the naturally varying position vectors of 
the satellites produce unique uvw sampling points, which consequentially improve the PSF. 

In order to avoid interference either the cluster must be located far from Earth-based REI and 
ionospheric distortions, such as Earth leading/trailing orbits and Sun-Earth Lagrangian points. However, 
by increasing the distance to the Earth, the communication with the Earth becomes more difficult. 
Alternatively, REI shielding can be achieved by positioning the array on the far side of the moon. The 
Radio Astronomy Explorer RAE-2 showed that interference at very low frequencies is reduced by 2 orders 
of magnitude behind the moon, making it an ideal location for radio astronomy observations (Alexander 
et ah, 1974). However, during the eclipse behind the moon the cluster has no communication with Earth. 
The following section discusses the quest to find an optimum balance between down-link data rates and 
maximizing observation time, emphasizing the challenges in various deployment locations. 

7.1 Lagrange points and Moon-farside 

The relative velocities of the satellites are minimal at Lagrangian points and hence these locations offer 
increased positional stability for longer time intervals. Therefore, the Lagrange points are an optimal 
choice to increase the integration time of the observations and also the mission lifetime. The Earth- 
Moon L4 and L5 are much closer for Earth based communication, however are suspected to be radio 
quiet relative to the Earth-Moon L2. The Earth-Moon L2 located at ~61347km away from the Moon, 
is still in the cone of radio-silence and is sufficiently shielded from REL However, this lagrange point 
may not be a favorable deployment location since transmission in this radio quiet zone may affect future 
missions (Maccone, 2005). The Sun-Earth L4 and L5 points are too far and subsequently limit downlink 
rates. In contrast, the Sun-Earth L2 liberation point at ~1.5 Million km away from Earth, is a tradeoff 
between downlink data rate, REI avoidance and increasing Tint- Although this is a stationary point, in 
practice a satellite operating at L2 will experience a gravity gradient with a slow and steady outward 


18 













Figure 12: Altitude vs Orbital period in Lunar orbit: A satellite cluster orbiting the moon over the far- 
side enters a “cone of silence” behind the moon once every orbit. During this phase, the moon shields the satellite cluster 
from Earth-based and solar interference, thus permitting relatively noise-free scientific observations. These eclipse periods 
can be extended by inserting the satellite network at a higher altitude, which also offers more orbital stability (Rajan 
et ah, 2011). However, a longer eclipse period behind the moon and subsequently larger science duty-cycle constrains the 
communications with Earth. 


drift. Such a scenario is preferred by the FIRST (Bergman et ah, 2009) and SURO-LC (Baan, 2012) 
studies. The SURO-LC proposes a array of 8 daughter satellites drifting slowly in Lissajous orbit and 
a mothership at a fixed distance of 10km from the cluster. While such a mission will provide enhanced 
imaging performance with improved uvw coverage and longer integration times, the downlink data rate 
is estimated to be 2 — 3 orders of magnitude less than an Moon based array using current technology 
(Rajan et ah, 2011). 

7.2 Orbiting the Moon 

An equatorial orbit around the Moon presents a relatively easier down-link to Earth and sufficiently long 
eclipse times behind the Moon w.r.t. Earth. The long eclipse time periods shield against radio noise from 
Earth and enable the science observations. In the DARIS study, to increase the predictability of the 
relative positioning, the reference orbit around the Moon was chosen to be circular which additionally 
also decreases the chance of collisions (Saks et ah, 2010). The array formation is build up from different 
relative orbits with a different relative inclination around a reference orbit (Eigure 12(a)), where each 
of these relative orbits contain several node spacecrafts (Eigure 12(b)). The reference orbit determines 
the duration of the eclipse time and subsequently the science duty cycle, which is increased by aligning 
an orbit with the Earth-Moon plane and/or by lowering the altitude (Boonstra et ah, 2011). As seen 
in Eigure 12(c), the Eclipse time period can be increased by decreasing the orbital altitude, however 
consequently the percentage of the orbit in the shade increases. In addition, by decreasing the orbital 
altitude, the relative range rates of the satellites also increase, which in turn affects the baseline stability. 
Hence, a balance between the relative velocity and the eclipse time must be found. When including the 
perturbations of the Earths’ gravity field, the irregularities in the lunar gravity field, the solar gravity 
field and the solar pressure, a constant drift of the relative orbits occurs. Coincidentally, this drift is 
mainly along the in-track direction of the reference orbit which can be compensated by adjusting the 
semi-major axis of the spacecraft node. In essence, a circular orbit in the Lunar equatorial plane offers 
a stable orbit, provided a trade-off is achieved between eclipse time and the satellite range rates. 

7.3 Orbiting the Sun 

A potential reference orbit for formation flying around the Sun is the Earth orbit itself. However, if 
the satellites are too close to Earth, then the terrestrial interference is a major disturbance to science 
observations. Alternatively, an orbit around the Sun with a different eccentricity than the Earth orbit 
keeps the satellite array at 4 to 10 million km from Earth, which is far enough to offer both stability and 
also reduce radio noise from Earth. The large distance separation severely limits the available down¬ 
link bandwidth upto at least an order magnitude compared to the Lunar orbits. In view of an optimal 
balance between increased data-downlink and REI free science observation, we choose the Earth orbit as 
a reference orbit with the satellite nodes orbiting at a distance of 4 to 10 million km from Earth. Hence, 
even though the constellation orbits the Sun as a central body, the reference orbit does go around the 
Earth, from leading to trailing. One of the many benefits of this particular orbit is that it is relatively 
stable for 10 years and allows continuous scientific observations. However unlike the lunar orbital design, 
this reference orbit is eccentric and highly sensitive to small changes due to large difference between the 
semi-major axis and the relative pairwise distances of the satellites (Boonstra et ah, 2011) (Saks et ah, 
2010). Eurthermore, the time period of the reference orbit is equal to the period of the relative orbits, 
which causes the formation to drift in all directions. With reference to the reference orbit, the cluster 
will slowly expand with time and hence offers unique sampling points for interferometry. One of the 
key advantages of this orbit is the low relative range rates which facilitates longer integration period. In 
addition, the relative range rates of the satellites in this orbit are less than 20cm/s for 3 years. Despite 
this advantage, the solar orbit is sensitive to small errors in velocity and the relative orbits are stable 
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only for change in iniection velocities upto O.lmm/s, which can be compensated using minor corrections 
(Saks et ah, 2010). 


8 Summary and Discussion 

A satellite cluster of less than 10 nodes is scientifically very interesting and meets the requirements 
for the extra-galactic survey science cases in terms of resolution and sensitivity. At least 4 antennas 
observing at 30 MHz for more than a year is sufficient to achieve the confusion limit of 65 mJy with 1' 
resolution, in which case over a million sources can be detected (Section 2.3). Moreover, even with fewer 
antennas, transient science cases such as bright Jupiter-like flares and Crab-like pulses can be addressed. 
All the satellites will be equipped with 2 (or 3) 5m dipole antennas (or two 2.5 monopoles) to observe 
the < 30MHz spectrum (Section 3). 

For a nominal observational bandwidth of > IMHz, each satellite is estimated to generate > 6 
Mbits/s, which must be correlated in space to minimize downlink data rate to Earth. In both centralized 
and distributed scenarios. The processing requirements for filtering and correlation is negligibly small for 
upto 50 satellites and can be readily incorporated into the On Board Computer (OBC) (Section 4). To 
establish the Inter-satellite link, satellites will be equipped with patch antennas to transmit the desired 
> 6Mbps data rate. The ISL budget analysis shows that in the Centralized scenario using 2.45GHz ISM 
band, the Node to mothership link can be established with 5W over 50km distance with a positive link 
margin. However, in the Distributed scenario upto 15W is desired to establish a link over 100km, which 
could be improved using clustering schemes and multi-hop communication (Section 5). Nonetheless, the 
proposed Distributed framework remains indispensable for large and scalable array of > 10 satellites, 
where SPOF must be avoided. 

The on-board clock on all satellites must have an Allan-deviation of < 10“^^, which can be met 
by a Rubidium or OCXO clock (Section 4.3). The current bulky space-qualified clocks, such as EADS 
OCXO, will potentially be replaced by light-weight and low-power on-chip atomic clocks e.g., SA.45s. In 
inaccessible (e.g., Moon-farside) or far-away deployment scenarios (e.g., Lagrange points), the satellites 
can be synchronized and localized using MGLS like algorithms, which enable the satellite network to 
be self-reliant co-operative network with minimum dependence on Earth-based ground stations (Section 
6). In addition, the orientation of the satellites can be estimated using the sensors in the Attitude and 
Orbit Gontrol System (AOGS) which include the sun sensor and star trackers. All satellites will also be 
equipped with sufficient propulsion to ensure precise deployment and to maintain the maximum baseline 
separation of 100km (Section 7). 

8.1 Technological challenges for ULW arrays 

The actual satellite implementation is intricately connected to the specific mission requirements, the 
number of satellites, the active choices in network architecture and the deployment location. However, 
recent studies which investigated centralized scenarios for an ULW array give insights into the current 
state-of-the-art space technology. Figurel3 shows the mass and power breakdown for the DARIS mission, 
where all subsystems use only existing and tested off the shelf components (Boonstra et ah, 2011). 
The power consumption for the daughter node and the mothership was estimated at 160W and 502W 
respectively. Reliable and highly efficient solar panels based on triple junction GaAs cells were employed 
on both the mothership and Daughter nodes to meet the power requirements. Furthermore, the dry 
mass of each daughter node was estimated at ~ 100kg and the mothership at ~550kg. 

In comparison to DARIS, futuristic missions such as OLFAR, are expected to be lighter by two 
orders of magnitude and consuming an order of magnitude lesser power (see Table 1). The reduced mass 
and power requirements will not only enable a larger array of antennas for radio astronomy, but can 
potentially enable the system to piggy-back on other missions, without the need for a dedicated launch 
vehicle. Thus, future missions will possibly consist of relatively cheaper nano-satellites with miniaturized 
and power-efficient subsystems. 

The intra-satellite communication between the satellite nodes is a fundamental bottleneck, which 
limits the bandwidth of observation and possibly the achievable baseline for radio astronomy imaging 
(see Section 5.2). In addition to limiting the feasibility of the science cases, the power consumption 
of existing technologies is also high. The DARIS project indicates > 25% of the power consumption 
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Daughter Node Mass Budget 



■ Link to Mothership (14kg) 

■ Harness (7kg) 

■ Payload (1kg) 

■ Structure/Mechanisms (18kg) 

■ AOCS(8kg) 

■ Propulsion (20kg) 

■ Power+ODBH (23kg) 

■ Solar Array (5kg) 

■ Thermal (4kg) 


Daughter Node Power Budget 


Mothership Mass Budget 



I Link to Daughters (11kg) 

I Link to Earth (49.5kg) 

I Harness (49.5kg) 

I Payload (5.5kg) 

■ Structure/Mechanisms (215kg) 

■ AOCS (27.5kg) 

I Propulsion (82.5kg) 

I Power (55kg) 
lODBH (22kg) 

■ Solar Array (16.5kg) 

Thermal (16.5kg) 

Mothership Power Budget 



■ Link to Mothership (45.2W) 

■ Harness (8.1W) 

■ Payload (21.1W) 

■ Power+OBDH (46.5W) 

■ Thermal (14.1W) 

■ AOCS (24.9W) 



■ Link to Earth (155.6W) 

■ Link to Daughters (60.2W) 

■ Harness (25.1W) 

■ Payload (85.3W) 

■ OBDH (35.IW) 

■ AOCS (60.2W) 

■ Thermal (40.2W) 

■ Power (40.2W) 


Figure 13: Mass and Power budget analysis of the DARIS mission: The DARIS mission consists 

of 8 Daughter nodes and 1 centralized mothership. The mass (and power) of each Daughter satellite and mothership was 
estimated to be 100kg (160W) and 550kg (502W) respectively. 


for communication for both the daughter node and mothership respectively (FigurelS). The satellite 
network to Earth communication limits the possible number of satellites in the cluster. Moreover, one of 
the limiting factors for the number of satellites is the downlink data rate of the satellite network, such as 
the HBA XDU for the centralized architecture (see Section 5.3). In case of distributed architecture, the 
satellite-swarm will employ diversity schemes to cooperatively downlink data to earth (Budianu et ah, 
2015). 

Further potential research areas identified during the study include the antenna design for observation 
frequencies below lOMHz, development of efficient imaging techniques for radio astronomy, high speed 
and robust RF Inter-satellite communications techniques (Budianu et ah, 2013) and investigating control 
and reliability of large satellite arrays (Engelen et ah, 2014). In addition, observability challenges such 
as the unknown REI environment at the desired deployment location must also be investigated, possibly 
by a > 2 satellite interferometer via a precursor mission. 

8.2 Conclusion 

The frequency window of < 30 MHz opens a new realm of interesting science cases and yet remains the 
last unexplored frequency regime in astronomy. To achieve the science objectives at these wavelengths 
with desired resolution and sensitivity, a dedicated space-based ULW array is necessary. Recent advances 
in technology and computing resources have improved both the feasibility and scientific desirability of 
such a space-based array. In this article, we justified the need for a space-based antenna array for ultra- 
long wavelength radio astronomy and discussed various subsystems needed to achieve the desired science 
cases. More recently concluded projects such as DARIS, EIRST, SURO have shown feasibility of such an 
array. In particular, the DARIS project showed that a cluster of less than 10 satellites can be launched 
using current off the shelf technology. An expanded set of science cases can be targeted by scaling the 
number of satellite nodes, extending the frequency range of observation and increasing the instantaneous 
bandwidth. However, this would significantly increase the mass, power consumption and eventually the 
cost of the mission. The on-going work on miniaturized nano-satellites may overcome this bottleneck 
and pave the way for feasible and affordable missions in the future. 
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